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ABSTRACT 


One objective of this thesis was to investigate the effect of details, such as the 
windows of high-rise buildings, on the radiowave propagation in the dense urban 
environment through modeling and simulations. If adding windows does not 
significantly change the signal distribution on average, it may not be necessary to build 
such a detailed model. Simulations are performed using several levels of detail and the 
results compared to estimate the impact of the fine details on the signal level. 

A second issue is base station antenna coverage. The antenna gain, half power 
beamwidth (HPBW), location, and pointing angle should be chosen to give the maximum 
coverage over a specified sector. Simulations can be used to select the optimum set of 
base station properties. Specifically this research looks at the coverage over a quadrant 
from the two sectored antennas versus a single one. 
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EXECUTIVE SUMMARY 


Since the proliferation of wireless devices in the 1990s, the study of radio 
frequency (RF) propagation in the urban environment has received much attention. With 
the inereased demand for products and services ineorporating wireless technology, the 
need for researeh into the behavior of RF propagation in the urban environment has come 
to the forefront. 

Cellular telephone service providers are installing base stations in urban 
environments at loeations that will ensure optimum customer aeeess. These loeations ean 
be determined by condueting a eostly RF survey of the coverage provided by the 
proposed base station location. Surveys may involve the plaeement of a transmitter and 
antenna at some proposed loeation, and then field strength is measured at multiple points 
throughout the serviee volume. This approach to the identifieation of base station 
plaeement is not only eostly, but time consuming as well. 

Computer simulation offers an alternative approaeh to the identification of 
optimum antenna placement. A computational electromagnetic (CEM) code, sueh as the 
Urbana Wireless Toolset, can be used to prediet the RF propagation in some specified 
urban region. The simulation approaeh is not only time-saving, but mueh less costly than 
the RF survey. 

One objeetive of this thesis was to investigate the effect of details, sueh as the 
windows of high-rise buildings, on the radiowave propagation in the dense urban 
environment through modeling and simulations. 

A seeond issue is base station antenna coverage. The antenna gain, half power 
beamwidth (HPBW), loeation, and pointing angle should be ehosen to give the maximum 
coverage over a specified sector. 

In this thesis, a group of city blocks roughly modeled after New York City 
Manhattan Island near Central Park was used for simulation. The buildings were adapted 
from an existing CAD model of the city downloaded from 3d Cadbrowser.com. Each 
building was then reshaped and resized to mateh the dimensions of the actual city blocks. 
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Based on computer simulations which include ray tracing augmented by edge 
diffraction, we conclude that adding windows does not significantly change the signal 
distribution, on average. Data analysis from simulations also indicated that varying 
frequency over the specified range would have little impact on the overall results. 

This thesis also analyzed the case where two antennas were used in combination 
and compared to the performance when only one larger beamwidth antenna was used 
with higher transmit power. The simulations results indicated that there was no apparent 
difference in performance between these two cases studied. 



1. INTRODUCTION 

A. BACKGROUND 

Since the proliferation of wireless deviees in the 1990s, the study of radio 
frequeney (RF) propagation in the urban environment has received much attention. With 
the inereased demand for products and services ineorporating wireless teehnology, the 
need for researeh into the behavior of RF propagation in the urban environment has eome 
to the forefront. 

Cellular telephone serviee providers are installing base stations in urban 
environments at loeations that will ensure optimum customer aeeess. These loeations can 
be determined by eonducting a eostly RF survey of the eoverage provided by the 
proposed base station loeation. Surveys may involve the plaeement of a transmitter and 
antenna at some proposed loeation, and then field strength is measured at multiple points 
throughout the serviee volume. This approaeh to the identifieation of base station 
plaeement is not only eostly, but time eonsuming as well. 

Computer simulation offers an alternative approach to the identifieation of 
optimum antenna plaeement. A eomputational eleetromagnetie (CEM) eode, sueh as the 
Urbana Wireless Toolset, ean be used to prediet the RF propagation in some speeified 
urban region. The simulation approaeh is not only time-saving, but mueh less eostly than 
the RF survey. 

A major eonsideration in seleeting the simulation approaeh is whether it is 
suffieiently aeeurate. The underlying eleetromagnetie (EM) theory is well known so it is 
a question of how aeeurately the EM equations are implemented and solved. Beeause of 
the frequeneies involved ( > 1 GHz ) high frequeney teehniques are usually employed. It 
has always been assumed that building wall refieetions and corner and roof edge 
diffraetions are the dominant propagation meehanisms. If this is the ease then very erude 
models of buildings would suffiee and would make the simulation approach much more 
appealing. However, reeent studies by Ghoraish [1] have shown that signal seattering 
from objeets sueh as lampposts, traffie lights, signboards ean be eomparable to the wall 
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reflections. Calculation of signals scattered from trees [2] have been done and the results 
were published recently. 

These seemingly minor details can have a significant effect on nearby observation 
points. However, these are not the only objects that contribute to the scattering of 
radiowaves in the urban environment. There are the fine building details: window 
frames, railings, stairs, etc. A recent paper [3] suggests that multiple diffractions from 
periodic grating structures like railing slats dominate at sizable distances from the 
structures. Further investigation and evaluation of different levels of detail in the urban 
models seems to be necessary. 

B. OBJECTIVES 

Over the years, there have been numerous studies of how radiowaves could be 
affected when propagating in an urban environment. Ray tracing is one of the methods 
that offer the potential for the greatest accuracy in complex environments and has been 
the focus of much research over recent years. Ray tracing includes the method of 
geometrical optics (GO) and the geometrical theory of diffraction (GTD). Ray tracing 
requires extremely accurate terrain and clutter databases and calculations can be 
computationally expensive due to the large number of rays that must be traced [4]. 

One objective of this thesis was to investigate the effect of details, such as the 
windows of high-rise buildings, on the radiowave propagation in the dense urban 
environment through modeling and simulations. If adding window and door frames does 
not significantly change the signal distribution on average, then it may not be necessary 
to build such a detailed model. Simulations are performed using several levels of detail 
and their results compared to estimate the impact of the fine details on the signal level. 

A second issue is base station antenna coverage. The antenna gain, half power 
beamwidth (HPBW), location, and pointing angle should be chosen to give the maximum 
coverage over a specified sector. Simulations can be used to select the optimum set of 
base station properties. Specifically this research looks at the coverage over a quadrant 
from two slightly overlapping sectored antennas versus a single one with a larger HPBW. 


2 



A test model was used to simulate the effeets of different material properties, 
antenna patterns, and frequeney variation. The simulations also aceounted for diffraetion, 
which is particularly significant in outdoor urban propagation environments. 

The simulation results predict the local mean power received at any given point 
on the observation grid. For each observation point, the vector sum of the multipath field 
can be computed. Data can then be presented graphically for quantitative evaluation or it 
also can be presented tabulated for more detailed examination. 

C. RELATED WORK 

Following are some recent studies that utilize the ray-tracing technique to model 
urban propagation. 

In [5] three diversity techniques were investigated to assess their performance in 
an urban environment. They were spatial diversity, polarization diversity, and angle 
diversity. These techniques were simulated for a UAV platform operating over a small 
city. The urban model of single building material was used with limited detail of the 
building. The results indicated that angle diversity has the most significant improvement 
in received signal strength as compared to the other two techniques. 

Pala [6] examined propagation for several frequencies. The results indicated that 
high frequencies suffer attenuation more rapidly in lossy buildings. Vertically polarized 
antennas generally give better results than horizontally polarized antennas. 

Ghoraishi [1] identified objects such as signboards, traffic signs, etc. that 
contributed to the scattering of the waves in the urban environment. This paper 
concluded that scattering from these objects can be comparable to wall reflection and, in 
any case, it has a significant contribution to the non-line of sight (LoS) received waves. 
Any prediction of the urban propagation channel needs to carefully consider and evaluate 
these scattering effects. Finally, [4] examined the effect of architectural features such as 
repeated rows of railings on propagation. It found that double diffraction can be the 
dominant mechanism in determining the field decay at relatively large distances (~ 0.5 
km) from the features when they are still in view. An example of such features is 
illustrated in Figure 1, which is a building in the French quarter in New Orleans. Similar 
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features appear on more modem buildings as shown in Figure 2, whieh is located 
adjacent to Central Park in New York City. 



Figure 1. Periodic stmcture of railings, French Quarter New Orleans (from [7]) 



Figure 2. Repeated rows of railings and balconies on an apartment building, NYC (from 

[ 8 ]) 


D. THESIS OUTLINE 


Chapter II presents concepts in electromagnetics that are relevant to the problem 
under consideration. They include propagation mechanisms and antenna fundamentals. 
Software tools such as the Urbana Wireless Toolset, and Rhinoceros Computer Aided 
Design (CAD) are also discussed in this section. A flowchart is presented to illustrate the 
steps involved in mnning the Urbana program. 
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Chapter III presents results from simulations and provides data analysis. Data is 
presented in both graphieal and tabular formats. Graphical data are used to demonstrate 
the effect of the details of the building on the radiowave propagation. The tabular data 
format is used for comparison between the different model conditions and variations in 
the base station antenna coverage. 

Chapter IV presents conclusions, suggestions for future research, and final 
remarks. 
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II. RADIO WAVE PROPAGATION IN THE URBAN 

ENVIRONMENT 

In this chapter, some eoncepts in eleetromagneties that are relevant to radiowave 
propagation in the urban environment are diseussed. Also diseussed are some of the 
eommon models used for the predietion of signal strength. The Urbana wireless toolset 
and Rhinoceros Computer Aided Design (CAD) software are also deseribed. 

A. PHYSICS OF PROPAGATION 

I. Electromagnetic Radiation 

Electromagnetic radiation is generally described as a propagating wave in space 
with electric and magnetic field components. These components oscillate at right angles 
to each other and to the direction of propagation, and they are in phase with each other in 
free space. 

Maxwell derived a wave form of the electric and magnetic equations, revealing 
the wavelike nature of electric and magnetic fields and their symmetry. According to 
these equations, a time-varying electric field generates a magnetic field and vice versa. 
Therefore, as an oscillating electric field generates an oscillating magnetic field, the 
magnetic field in turn generates an oscillating electric field, and so on. These oscillating 
fields together form an electromagnetic wave [9]. 

Modern electromagnetism is based on a set of four fundamental relations known 
as Maxwell’s equations. Using phasor quantities with a e^“"time convention. Maxwell’s 
equations have the form 


vb = p^ 

2-1 

VxE = -jcoB 

2-2 

VB = 0 

2-3 

V X // = J + jcoD 

2-4 
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where E and D are the electrie field intensity and eleetric fiux density, respeetively. 
They are related hy D = sE , with s being the permittivity of the material. The magnetic 
field intensity, H , and magnetic fiux density, B , are interrelated hy B = /uH , with // 
being the magnetic permeability of the material. The charge density per unit volume is 
and the volume current density is J [9]. 

Maxwell’s equations for a charge free medium can be used to obtain the wave 
equation for the electric field 

V^E-Y^E = h 2-5 

where y = co-sf/us is the propagation constant andct> = In f. 

Two important solutions to the wave equation are plane waves and spherical 
waves. A plane wave has the general form 

E{R) = 2 - 6 

where R = hc + yy + zz is a position vector to the point (x,y,z). The vector E^ 
determines the polarization of the plane wave. The vector propagation constant is 


r = r.x+ryy+rj 2-1 

where the components are given by 

= coyfjus y-x 2-8 

ry=o)^Y-y 2-9 

Y^ = co^fpe y-z, 2-10 


and y is a unit vector in the direction of the propagation. 

Spherical waves originating at the origin have the form 

-yr 

E{r) = {E,e + E^h^. 2-11 

r 
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Both plane and spherical waves are transverse electromagnetic (TEM) waves. The 
magnetic field intensity is obtained from 


H{r) = 


/xE(r) 

ri 


2-12 


The impedance of the medium is = 

The properties of an electromagnetic wave, such as its phase velocity, , and 

wavelength, X , are governed by angular frequency, co , and the three constitutive 
parameters of the medium; s, /u , and cr, such that 



u 


p 


1 

4ms’ 


2-13 



2-14 




s = s^s^=s^{s’r-jsl) = s^ 


^r-J- 


(J 


cos, 


0 y 


2-15 

2-16 


If the medium is lossless (cr = 0 and ju” = 0 )> the wave does not suffer any 
attenuation as it travels through the medium. However, when the wave is propagating 
through a lossy medium, (cr > 0 or 0), absorption occurs as the wave passes 
through. Depending on the material the losses can be due to induced polarization, 
magnetization, or ohmic loss. 

2. Relationship Between Spherical Waves and Plane Waves 

When energy is emitted by a source, such as an antenna, it expands outwardly 
from the source in the form of spherical waves, as depicted in Figure 3 for an observer in 
the far field. The spherical wave travels at the same speed in all directions and expands 
at the same rate. If an observer is located sufficiently far from the source (R^ co), the 
wavefront of the spherical wave appears approximately planar, as if it were part of a 
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uniform plane wave with uniform properties at all points in the plane tangent to the 
wavefront as shown in Figure 3 [9]. Therefore, even though pure plane waves do not 
exist, the plane wave approximation can be used to accurately model real-world 
problems. 



Figure 3. Spherical wave (from [11]) 

B. RAY OPTICS 

Geometric optics (GO), also known as ray optics, represents an electromagnetic 
wave by a ray denoting the direction of travel of the wavefront. Neither the phase nor the 
polarization of the wave is accounted for explicitly in geometric optics. Thus, geometric 
optics is a ray approximation of physical optics [9]. 

1. Geometric Optics 

In general, the GO method provides fairly accurate results whenever the size of 
the object illuminated by an electromagnetic wave is much larger than the wavelength, 
/I. GO also describes wave behavior upon reflection or refraction at an interface 
between two materials, based on the following assumptions [11]; 

1. Wavefronts are locally plane and waves are TEM. 

2. The wave direction is specified by the normal to the equiphase surface. 
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3. Rays travel in straight lines in a homogeneous medium. 

4. Polarization is constant along a ray in an isotropic medium. 

5. Power in flux tube is conserved as depicted in Figure 4. 

jj Jj 2-17 

Areal Areal 



Figure 4. Flux tube (from [11]) 


6. The reflection of rays obeys Snell’s law. 

7. The reflected field is linearly related to the incident field at the reflected 
point by a reflection coefficient. 

2. Geometric Theory of Diffraction 

The geometric theory of diffraction (GTD) is a supplement to GO by the inclusion 
of diffraction. GTD eliminates some of the problems of GO, namely that GO does not 
account for the edge diffracted field. The total field at an observation point P can be 
found by the sum of GO and diffracted components [11]; 

E{P) = %oiP) + E,{P). 2-18 

Figure 5 illustrates how edge diffraction contributes to the field at an observation 
point P. The reflected and refracted rays comprise the GO part of the field. Multiple 
reflections are also possible although none are shown. Only single diffractions are 
shown; that is SAP and SBP. There are also multiple diffractions such as SABP. Finally, 
there are mixed terms such as reflected - diffracted that are less significant than the GO 
and GTD terms, and can be ignored. 
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Figure 5. Principle of ray tracing model 


C. ANTENNA FUNDAMENTALS 

Antennas are treated as either transmitting or receiving as appropriate for the 
particular situation. In the receiving mode, antennas act to collect incoming waves and 
direct them to a common feed point where a transmission line is attached. Antennas have 
directional characteristics; that is, electromagnetic power density is radiated from a 
transmitting antenna with intensity that varies with angle around the antenna. 

1. Antenna Radiation Pattern 

The radiation pattern of an antenna gives the angular variation of radiation or 
reception at a fixed far field distance from the antenna. The antenna responds to an 
incoming wave from a given direction according to the pattern value in that direction. 
The normalized antenna pattern factor, F{6, , describes the shape of the radiation 

pattern in three-dimensional space. For a short dipole oriented along the z axis the 
normalized antenna pattern is 

F{6,(l)) = ?,m6. 2-19 

The complete three dimensional radiation pattern for the ideal dipole resembles a 
“doughnut” pattern as shown in Figure 6. The £'-plane and //-plane patterns of the ideal 
dipole are shown in Figure 7 [9]. For a vertical dipole as shown, the pattern cut 
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transverse to the dipole axis is referred to as the azimuth pattern. The pattern cut that 
contains the dipole (z-axis) is referred to as the elevation pattern. 



Figure 7. (a) Ideal vertical dipole, (b) H-plane (azimuth) radiation pattern, (c) E-plane 

(elevation) radiation (from [9]). 


2. Directivity and Gain 


The directivity of an antenna is defined as the ratio of the radiation intensity in a 
certain direction to the average radiation intensity, or 


D = 


1 


4;r 


2-20 


where is referred to as the beam solid angle. 
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For an antenna with a single main lobe pointing in the z-direction as shown in 
Figure 8, Q^may be approximated as the product of the azimuth and elevation half 

power beamwidths ^„pBw [9] 



^HPBW^HPBW 
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(b) Half-power beamwidth 
m x-z plane 



(c) Half-power beamwidth 
m ,v-r plane 


Figure 8. Antenna directivities (from [9]). 

The gain, G, of an antenna is defined relative to the power into the antenna and 
thus includes loss. In practice, the gain of an antenna can be determined from its effective 
area, , and the wavelength, X, by 


G = ^D = 




2-22 


where ^ is the radiation efficiency: 




radiated 


input 


( 0 <^< 1 ). 
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In this thesis, the antenna radiation is modeled as a cos'" 6 power pattern. The 
field pattern is of the form 
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E{6,(j)) - -(6* cos"" ^ cos cos'" 6’sin^zJ), 

r 



Following is the sample caleulation for the ERP. If the HPBW for eaeh antenna is 
given then its corresponding gain can be eomputed from the speeified HPBW. The 
HPBW is twice the half power angle 


^HPBw 2<!eos 


( 

0.5'” 

V 
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From the HPBW it is possible to solve for m 

log(0.5) 


m = 


r 1 

C 0 


1 

eos 

1 ^HPWB 


1 

1 ' 

1 2 

7 

1 
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For =55 deg , m = 5.78 and = 110deg , m = 1.247 


The direetivity is [10] 

D = 2{m + \) 2-27 

and the ERP for a lossless antenna (G = Z) ) is 

ERP = P,G. 2-28 


3. Antenna Polarization 

The energy radiated by any antenna is eontained in a transverse eleetromagnetie 
wave that is eomprised of eleetrie and magnetie fields. Eor plane and spherieal waves 
these fields are always orthogonal to one another and orthogonal to the direetion of 
propagation. The eleetrie field of the eleetromagnetie wave is used to deseribe its 
polarization and henee, the polarization of the antenna. In general, all polarizations fall 
under the general ease of elliptieally polarized. In this general ease, the total eleetrie field 
of the wave is eomprised of two linear eomponents, whieh are orthogonal to one another. 
Eaeh of these eomponents has a different magnitude and phase [9]. At any fixed point 
along the direetion of propagation, the total eleetrie field would traee out an ellipse as a 
funetion of time. 
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The two limiting cases of elliptical polarization are circular and linear. A 
circularly polarized electromagnetic wave is comprised of two linearly polarized electric 
field components that are orthogonal, have equal amplitude, and are 90 degrees out of 
phase. In this case, the polarization ellipse traced by the wave is a circle. Depending 
upon the direction of rotation of the electric field vector, the wave will be left-hand 
circularly polarized or right-hand circularly polarized. The phase relationship between 
the two orthogonal components, +90 degrees or -90 degrees, determines the direction of 
rotation. A linearly polarized electromagnetic wave is comprised of a single electric field 
component, and the polarization ellipse traced by the wave is a straight line [9]. Most 
commercial wireless systems are designed to operate with linear polarization, usually 
vertical with respect to the ground. 

4. Free-Space Link Equation 

The two antennas shown in Figure 9 are part of a free-space communication link, 
with the separation between the antennas, R, being large enough for each antenna to be in 
the far-field region of the other [13]. The transmitting and receiving antennas have 
effective areas At and Ay and radiation efficiencies, and respectively. The antenna 
gains are related to their effective areas by Equation 2-22. The received signal power in 
free space can be expressed by the Friis transmission equation as follows: 


R = 


{AnRf 


2-29 


Free space path loss is the spreading loss in signal between two isotropic antennas 
(= 1), and it can be expressed as: 


L 


'fs 


Py _ 

P, (ATTRf ■ 
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In urban propagation applications it is common to use d in place of R for distance from 
the transmitter. 
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Figure 9. Transmitter-receiver configuration. 


D. THEORETICAL MODELS FOR URBAN PROPAGATION 

The urban and suburban propagation problems are complicated because the fields 
in the immediate vicinity of the portable or mobile radio are a superposition of localized 
multipath scattering. The signal strength varies from peak levels of a few decibels above 
the mean or median level to tens of decibels below the peaks in deep fades. 
Consequently, we often rely on a statistical description of the signal levels in the vicinity 
of the observer that states the local average and a description of the variation. There are 
many excellent models available for evaluation, however, only two models are 
considered; the diffracting screens model and the COST 231 model. 

I. The Diffracting Screens Model 

The urban propagation environment is complex and varied from one city to 
another, and even within a city. Simple approximations for the environment that are 
computationally efficient and accurate for a wide range of problems are desired. 

The Walfisch-Bertoni [14] model is used widely for predicting the average path 
loss for mobile systems in urban areas. The model assumes that the street grid in a typical 
city organizes buildings into rows that are nearly parallel and that an idealized 
representation for the urban environment would be as shown in Figure 10, where the 
precise height and spacing of the buildings have been ignored and the profile is 
characterized by just two parameters: the mean building spacing s, and the mean building 
height b. The transmitting antenna is positioned at height Hb above ground level and the 
receiver antenna is positioned at Hm relative to the ground [13]. With this simple 
representation, the average path loss from the transmitter to the receiver is the sum of free 


17 







space wavefront spreading, multiple forward-diffraction past the rows of buildings and 
diffraction over the final rooftop down to the receiver [14], 



Figure 10. Ray diffraction roofiine to street level (from [13]) 

The average excess loss above free space path loss, predicted by the 

Walfisch-Bertoni model can be written as the summation of two independent terms: the 
multiple diffraction loss as the field propagates past the rows of buildings, , and the 

diffraction loss from the last rooftop to the receiver, 

L =L,+L,. 2-31 

ave-ex md rtr 

2. The COST 231 Model 

The European Research Committee COST 231 created an urban propagation 
model based on the work of Walfisch and Bertoni and Ikegami [13]. The basic COST 231 
model uses Walfisch-Bertoni results to calculate urban environment propagation 
prediction along with Ikegami’s correction functions for taking into account the street 
orientation. The model was applied to the 800 MFlz to 1.8 GHz bands and tested in the 
German cities Mannheim and Darmstadt. The street orientation influence was found to be 
minimal while it was found to require considerable improvements when the antenna was 
at or below rooftop level. The accuracy of the COST-231 model is best for antenna 
heights much greater than the building heights, and deteriorates when the height of the 
antenna is less than the height of the buildings [15]. The model becomes inaccurate when 
the terrain topography is non-fiat or if the coating material of the land is inhomogeneous. 
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E. EMPIRICAL MODELS FOR URBAN PROPAGATION 


Many empirical models have been developed over the years to overcome some of 
the limitations and assumptions presented by the theoretical models. The empirical 
models are based on measured data from which curve-fitted equations are obtained to 
model propagation in areas of definable urbanization. Often, empirical models are city- 
specific and are tied to urban land use maps. Following are examples of a few empirical 
models that are available [13]. 

1. The Okumura Signal Prediction Method 

The experiment that Okumura [13] carried out was to measure the signal strengths 
in the vicinity of the city, Tokyo, over a wide range of frequencies, several fixed-site 
antenna heights, several mobile antenna heights, and over various irregular terrains and 
environmental clutter conditions. A set of curves relating field strength versus distance 
for a range of fixed-site heights at several frequencies was generated. Various behaviors 
in several environments, including the distance dependence of field strength in urban 
areas, and urban versus suburban differences were then extracted from these curves [13]. 

2. The Hata and Modified Hata Formulas 

The Hata empirical model uses a propagation equation split into two terms: a term 
that has a logarithmic dependence on distance and a term that is independent of distance. 
The Hata model also includes adjustments to the basic equation to account for urban, 
suburban, and open area propagation losses [16]. The Hata equation for propagation loss 
in an urban area is given by: 

=69.55 + 26.161og,o(/)-13.821og,o(AJ+ ^ 

[44.9 - 6.55 logio (/?,)] logio {d) + a^{hj 

where / is frequency in MHz, d is distance in kilometers, and \ is the base station 
height in meters. A mobile height correction function, {h ^), is applied for mobile 
antenna heights. In a medium city, Hata’s mobile height correction takes the form 

= [0.7-1.121og,o(/)]A„ +1.561og,o(/)-0.8 2 - 33 
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while in a large city, at 200 MHz and below, 


a,(AJ = l.l-8.291og?o[1.54AJ 


2-34 


while in a large city, at 400 MHz and below. 


a4(/^J = 4.97-3.21ogfo[11.75AJ. 


2-35 


3. Propagation Near Buildings 

Radiowave propagating in the urban environment can encounter additional losses 
due to the interference from nearby cells where the same frequency is reused. For short 
distances, there is a breakpoint beyond which the propagation law has increased losses 
with distance [13]. The following equation 2-36 can be used to explain the breakpoint 


2-ray 


Sin 


d 


2kd 


{kXdf 


2-36 


where for vertical polarization X is 

for horizontal polarization X is 

X = ^£^-cos\0) 

and; 

0 is the incident angle, 

£^ is the dielectric constant of the ground. 

Hi is the transmitter’s antenna height, 

H 2 is the receiver’s antenna height, 
d is the distance between antennas. 


2-37 


2-38 
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k is the wave number. 


At short distanees, power falls off proportional to . However, for longer distanees, the 
signal strength is falling off at an rate. The breakpoint oeeurs at the rate’s transition 
points from d^ to d"^. Figure 11 demonstrates the two-ray path model where //i = 4 m, 
H 2 = 1.6 m, and frequeney =1.9 GHz. The breakpoint oeeurs at a distanee of about 160 
m for both the vertieal and horizontal polarization [13]. 


dB. two ra> path 



Figure 11. Two-ray path propagation model showing the breaking point (from [13]) 

F. SOFTWARE TOOLS 

1. Rhinoceros CAD Software 

Rhinoceros CAD for Windows platforms is a software paekage that ean be used 
to draw and manipulate the geometry. For eventual use by Urban the file must be saved 
in 3-dimensional studio format (*.3ds). Sinee Urbana will read only a faeet file format, 
the 3ds format files must first be eonverted to faeet file format through the use of 
Menelaus tools. The next step is to extraet edges from the faeet file utilizing Menelaus 


21 




tools. Figure 12 shows the top level proeedure in flow ehart format to obtain the 
geometry input file, preparing an Urbana input file and post-proeessing using lurbana. 



Figure 12. Urbana Wireless Toolset and auxiliary programs (from [12]) 

The Urbana GUI ean present these results in graphie format or it ean also 
export for plotting and analysis using Matlab. One way to determine if there is any 
impaet due to windows is to ealeulate the differenee in the eleetrie field between a plain 
baseline model and detailed test model. Matlab software ean also be used to determine 
the average, standard deviation, and maximum field strength at eaeh point in the 
observation grid. These statisties are not normally provided by Urbana. 

2. Urbana Wireless Toolset 

The Urbana Wireless Toolkit is a eomputational eleetromagnetie tool for 
simulating wireless propagation in eomplex environments. The ray-traeing engine of the 
toolset is eoupled with proprietary algorithms to implement physieal opties, geometrieal 
opties, and diffraetion physies in produeing a three-dimensional (3-D) simulation. 

The ray-traeing proeess prediets the loeal mean power reeeived at a given point at 
X, y, and z. For eaeh point the veetor sum of multipath power is eomputed. The 
ealeulations also inelude the effeets of frequeney, polarization, material properties, and 
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varying antenna patterns. If desired, the model can include the effects of ray diffraction. 
The following list is the key input parameters for Urbana simulation software: 

1. Building facet and edge fdes. 

2. User defined building materials. 

3. Antenna patterns, placement, strength, frequency, and polarization. 

4. Observation grid. 

This toolkit contains two primary parts: Urbana_comp and lurbana. 
Urbana_comp is the stand-alone back-end application to perform the simulation and 
post-processing. In addition to executing the simulation jobs, its primary value is that it 
allows batch job processing without the need for a GUI interface to launch each job. 
There are two additional applications, Urbana and Urbana_rp, which are controlled by 
Urbana_comp and are transparent to the user. These applications run the analysis 
calculations when provided with data from Urbana_comp [17]. 

G. SIMULATION MODEL 

In this thesis, a group of city blocks roughly modeled after New York City 
Manhattan Island near Central Park was used for simulation. The buildings were adapted 
from an existing CAD model of the city downloaded from 3d Cadbrowser.com. Each 
building was then reshaped and resized to match the dimensions of the actual city blocks. 
The city model occupies an area of about 600 m by 500 m. The tallest building is at a 
height of 118 m. There are 24 separate city blocks that made up the simulation area. 

In order to determine if windows would have any significant impact on the 
radiowave propagation, two models were built for comparison. The two models were 
identical in size, number of buildings, and building materials. Only one of these models 
has windows on two blocks of buildings and is referred to as the detailed or test model. 
The other model is used as reference and referred to as the baseline model. Figures 13 
and 14 show the baseline model and test model respectively. Windows have been added 
to two blocks of buildings. The number was limited by the amount of computer memory. 
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Figure 13. Baseline model (a) eomplete view (b) detailed view 



Figure 14. Test model (a) complete view (b) detailed view 


This chapter discussed theory and concepts relevant to the study of RF 
propagation in the urban environment. The next chapter presents the results of some 
simulations of RF propagation in a modeled urban environment. A few city blocks of 
Manhattan Island in NYC were modeled to study the RF propagation in this environment. 
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III. URBANA SIMULATION RESULTS 


The previous ehapter diseussed coneepts that were relevant to the study of 
radiowave propagation in the urban environment. This ehapter will present the results of 
eomputer simulations that were eondueted to satisfy two objeetives: 

Objeetive 1 is to investigate the effeet of details, sueh as the windows of high-rise 
buildings, on the radiowave propagation. 

Objeetive 2 is to illustrate how the base station antenna pattern affeets eoverage. 

This ehapter begins with an explanation of the parameters used in the set up of 
simulations and follows with presentation of the data generated from these simulations. 

A. SIMULATION MODELS AND PARAMETERS 

A group of eity bloeks modeled after Manhattan Island in NYC was used for this 
simulation. These models were previously diseussed in Chapter II. NYC was ehosen 
beeause it is a good example of a dense urban environment. Table 1 shows the 
transmitting eharaeteristies of the antennas. Table 2 shows the input parameters for the 
simulations. The observation grids were set to an area of interest where buildings with 
windows are loeated. An observation grid is a region within the NYC model and it be 
defined in Urbana prior to the simulation. The height of the observation area and the 
resolution of eaeh eell also have to be predefined. Figure 15 depiets the observation grid 
for this simulation, and it is highlighted in blue. Figure 16 shows the loeation of the 
antenna. Figure 17 shows a elose up view of windows. 

The transmit antenna patterns and the eorresponding eoverage areas are shown in 

Figures 18, 19, and 20. The antenna loeation is at the top of the tallest building at the 

eomer of the model. Antennas 1 and 2 together provide eomplete eoverage; antenna 3 by 

itself provides eomplete eoverage. The z axis of the antenna (see Figure 8) is pointing 

downward slightly to give approximately uniform eoverage. Before the antenna is tilted 

the polarization is vertieal (z-direeted eleetrie field in Figure 8). Beeause the tilt angle is 

small the resulting polarization is very nearly vertieal. The antenna pattern surfaee mesh 

in the figures is sized for viewing, and the interseetion with the ground does not 
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necessarily represent the coverage area. The calculation of the antenna parameters is 
discussed in Section C. 

Table 1. Transmitting antennas characteristics 


Antenna 

HPBW 

(deg) 

Transmitting Power 
(W) 

Anteima Gain 
(dBi) 

1 

55 

1 

13.3 

2 

55 

1 

13.3 

3 

no 

3 

4.8 


Table 2. Simulation input parameters 


Input Parameter 

Input Values 

Observation height 

2 m 

Grid cell resolution 

0.5 X 0.5 m 

Frequency 1 

850 MHz 

Frequency 2 

900 MHz 

GO Ray-bounces 

5 

Diffraction 

single 



Figure 15. Observation grid (blue area) in over head view 
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Typical cell phone parameters were used for these simulations. Frequeneies of 
850 MHz and 900 MHz were ehosen. These frequeneies are 50 MHz apart, whieh is a 
typieal bandwidth for wireless modulations. 

In order to investigate the effeet of materials on propagation, two separate models 
were built. One of the models was built with all PEC material, and the other model was 
built with a more realistic building material, sueh as conerete. Instead of eutting holes in 
buildings for windows, radar absorbing material (RAM) was used behind the glass 
windows to simulate the effeet of radiowaves propagating through the windows and thus 
not refleeted baek out. This approximation eliminates the need to go baek to Rhino and 
out window openings in all of the buildings, whieh is a time oonsuming task. The 
material of the ground plane is the same as for all buildings in the model. The same 
thickness value is used for all building walls and ground plane. Table 3 summarizes the 
materials and their properties. 


Table 3. Properties of the Building materials 


Material 

Permeability 

Permittivity 

Thickness 

t 

A 

ft 

A 

s' 


(cm) 

PEC 

1 

0 

1 

0 

300 

CONCRETE 

10.2 

0.5 

1 

0 

300 

GLASS 

3.8 

0 

1 

0 

0.35 

RAM 

10 

100 

10 

100 

10 
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Figure 16. Complete model and observation grid three dimensional view 



Figure 17. Close up view of Windows 
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Figure 18. Antenna 1 pattern 



Figure 19. Antenna 2 pattern 
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Figure 20. Antenna 3 pattern 


B. SIMULATION RESULTS FOR BUILDING MATERIALS 

The goal of the first objective was to investigate the effect of details, such as the 
windows of high-rise buildings, on the radiowave propagation in the dense urban 
environment. In order to satisfy this objective, a set of reference data must be 
established. This is accomplished by performing a simulation with the baseline model 
that lacks the details. The simulation was then repeated with a more detailed test model. 
These simulations were performed with identical set up for all parameters. After the data 
was collected, an analysis was performed to determine the difference in signal levels 
between these models. Table 4 presents the results from the simulations where materials 
and test conditions were kept constant. Only the reflected and diffracted fields were 
considered. 
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Table 4. Simulation results of the diffraeted and refleeted fields for eonerete buildings 


Building 

Test 

Freq. 

(MHz) 

Ant 

Sample 

Baseline model 
(dBm) 

Detailed model 
(dBm) 

Delta 

Material 

Size 

Avg. 

Std. 

Max. 

Avg. 

Std. 

Max. 

H 


1 

23410 

-56.61 

2.30 

-42.67 

-55.91 

1.59 

-41.36 

0.71 

850 

2 

23410 

-55.86 

3.16 

-41.07 

-57.61 

4.92 

-41.40 

-1.75 

Pi 


3 

23410 

-55.26 

2.34 

-40.44 

-56.74 

3.82 

-41.86 

-1.48 

u 

z 


1 

23410 

-57.02 

1.89 

-44.00 

-56.23 

1.09 

-41.47 

0.79 

o 

u 

900 

2 

23410 

-55.92 

3.23 

-41.82 

-56.84 

4.15 

-41.57 

-0.92 



3 

23410 

-55.42 

1.69 

-44.56 

-55.86 

2.14 

-40.54 

-0.45 


The eolumn “Avg.” is the average power over all eells in the observation grid for 
an isotropie reeeive antenna. The overall average value of the baseline model is 
eonsistent with the range of the values in detailed model. These values are well above 
the typieal reeeiver’s threshold deteetion (-90 dBm), therefore this differenee will not 
have a signifieant impaet on the reeeiver’s signal at the observation area. 

The building material was ehanged from eonerete to simple perfeet eleetrieally 
eondueting (PEC) material, and the simulations were repeated. The results for PEC 
buildings are shown in Table 5. The all PEC buildings model is not realistie, but it 
served as a referenee for eomparison. The overall average is about 10 dBm higher for 
buildings with PEC material than the buildings with eonerete material. 

The magnitude differenee between the baseline and detailed models shown in the 
eolumn labeled “Delta” in Table 5 is small. This implies that the number of windows 
that were used in these simulations did not, on average, have signifieant effeet on 
radiowave propagation in the observation area. These results also show that ehanging 
frequeney from 850 MHz to 900 MHz did not have a signifieant effeet for either material. 
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Table 5. Simulation results of the diffraeted and refleeted fields for PEC buildings 


Building 

Material 

Ereq. 

(MHz) 

Ant 

Sample 

Size 

Baseline model 
(dBm) 

Detailed model 
(dBm) 

Delta 

Avg. 

Std. 

Max. 

Avg. 

Std. 

Max. 

PEC 

METAL 

850 

1 

23410 

-48.01 

2.38 

-33.74 

-48.49 

2.85 

-33.79 

-0.48 

2 

23410 

-47.17 

2.57 

-32.51 

-46.65 

2.05 

-32.47 

0.53 

3 

23410 

-45.20 

1.94 

-32.97 

-45.60 

2.34 

-32.62 

-0.40 

900 

1 

23410 

-48.63 

2.45 

-35.79 

-48.76 

2.59 

-35.47 

-0.13 

2 

23410 

-47.72 

2.53 

-34.30 

-47.34 

2.15 

-34.32 

0.39 

3 

23410 

-45.68 

3.15 

-32.81 

-45.07 

2.54 

-31.29 

0.61 


Figures 21 through 32 are the graphieal representations of the delta signal 
eontours. Eaeh figure presents a simulation eondition, and these eonditions are 
summarized in Table 6. The upper and lower limits of values displayed in eaeh figure are 
given in the eolor bars by the numbers on the right hand side. 

Table 6. Summarized simulation parameters 


Eigure number 

Antenna 

Building material 

Frequency 

(MHz) 

Figure 21 

1 

Concrete 

850 

Figure 22 

2 

Concrete 

850 

Figure 23 

3 

Concrete 

850 

Figure 24 

1 

Concrete 

850 

Figure 25 

2 

Concrete 

850 

Figure 26 

3 

Concrete 

850 

Figure 27 

1 

PEC 

900 

Figure 28 

2 

PEC 

900 

Figure 29 

3 

PEC 

900 

Figure 30 

1 

PEC 

900 

Figure 31 

2 

PEC 

900 

Figure 32 

3 

PEC 

900 


32 









Figure 21. Signal difference between baseline and detailed models for antenna 1 and 

concrete buildings/= 850 MHz 



Figure 22. Signal difference between baseline and detailed models for antenna 2 and 

concrete buildings/= 850 MHz 



Figure 23. Signal difference between baseline and detailed models for antenna 3 and 

concrete buildings/= 850 MHz 
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Figure 24. Signal difference between baseline and detailed models for antenna 1 and 

concrete buildings/= 900 MFlz 
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Figure 25. Signal difference between baseline and detailed models for antenna 2 and 

concrete buildings/= 900 MFlz 
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Figure 26. Signal Difference between baseline and detailed models for antenna 3 and 

concrete buildings/= 900 MHz 
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Figure 27. Signal difference between baseline and detailed models for antenna 1 and 

PEC buildings/= 850 MHz 
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Figure 28. Signal difference between baseline and detailed models for antenna 2 and 

PEC buildings/= 850 MHz 
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Figure 29. Signal difference between baseline and detailed models for antenna 3 and 

PEC buildings/= 850 MHz 
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Figure 30. Signal difference between baseline and detailed models for antenna 1 and 

PEC buildings/= 900 MHz 



Figure 31. Signal difference between baseline and detailed models for antenna 2 and 

PEC buildings/= 900 MHz 



Figure 32. Signal difference between baseline and detailed models for antenna 3 and 

PEC buildings/= 900 MHz 
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Figures 33 and 34 show the signal contours of a simulation for the entire model 
area. These figures also illustrate some effect of the individual components of ray 
tracing. A vertical dipole was used as the transmit antenna. 



(a) (b) 


Figure 33. Simulation of an entire area with antenna 3 (a) diffracted fields (b) reflected 

fields 


m 


dB V/m 



Figure 34. Simulation of an entire area for total field with antenna 3 (incident, reflected, 

and diffracted fields) 

Results from these simulations show that it is not necessary to build a model to 
the level of detail such as adding windows. Even though the signal contours show some 
variation of reflections around the area of the windows, it is not high enough to be 
considered significant. Therefore details of modeling should be focus somewhere else, 
such as billboards, and street signs. 
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C. SIMULATION RESULTS FOR ANTENNA PATTERNS 

Next, we examine how the base station properties (the antenna gain, half power 
beamwidth, loeation, and pointing angle) affeet the eoverage. A series of simulations 
was eondueted in support of this objeetive. Table 7 eontains results from these 
simulations. In this ease the total fields were used, beeause they are the aetual signals 
available to a mobile user. 

Table 7. Simulation results of the total fields for detailed model 


Material 

Freq. 

(MHz) 

Sample 

Size 

Combine Ant 1 & 2 
(dBm) 

Ant 3 
(dBm) 

Delta 

Avg. 

Std. 

Max. 

Avg. 

Std. 

Max. 

CONCRETE 

850 

23410 

-53.22 

2.66 

-38.22 

-56.74 

3.63 

-41.86 

3.52 

900 

23410 

-53.34 

3.02 

-38.56 

-55.86 

4.27 

-40.54 

2.52 

PEC 

850 

23410 

-44.25 

2.86 

-29.49 

-45.6 

2.14 

-32.62 

1.35 

900 

23410 

-44.67 

2.32 

-30.72 

-45.07 

3.07 

-31.29 

0.4 


In this simulation, antenna 1 and antenna 2 transmitted separately, and their fields 
were eombined to obtain the sum of the signal strength at eaeh point in the observation 
grid. The proeedure used was first to establish an effeetive radiated power (ERP) that is 
equal for both antennas under eonsideration. Therefore it was neeessary to eompute a 
transmitter power. Pi, for eaeh antenna used in these simulations that would produee the 
desired ERP. See Chapter II Seetion C for details. 

Eor antennas 1 and 2 the transmit power is 1 W (30 dBm) into eaeh antenna. A 
similar ealeulation was done to eompute Pt for antenna 3 that would yield the same ERP, 
whieh was found to 34.8 dBm. 

The data in Table 7 is based on both antennas 1 and 2 transmitting simultaneously 
eaeh with the full transmitting power. The average power of the eombined antennas 1 
and 2 is greater than the average value of the antenna 3, as expeeted. 

The results from these simulations demonstrated that it is not neeessary to foeus 
on the details of modeling and simulation, or the number of antennas needed for the 
simulation proeess. The average values over the observation grid are relatively 
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insensitive to these things. Perhaps attention should foeus more on billboards, signal 
signs, parked ears, or some other objeets that eould eause higher refleetion. 

The following ehapter presents the eonelusions for this thesis and offers a few 
topies for possible future work. 
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IV. CONCLUSIONS AND FUTURE WORK 


A. SUMMARY 

One objective of this thesis was to investigate the effect of details, such as the 
windows of high-rise buildings, on the radiowave propagation in the dense urban 
environment. A second objective was to investigate base station antenna coverage, that 
is, the antenna gain, half power beamwidth (HPBW), location, and pointing angle to give 
the maximum coverage over a specified sector. 

The buildings were adapted from an existing CAD model of NYC, Manhattan 
Island. Each building was then reshaped and resized to match the dimensions of the 
actual city blocks. Two models were built for comparison. The two models were 
identical in size, number of buildings, and building materials. Only one of these models 
has windows. The other model is used as a baseline for comparison. 

Urbana Wireless Toolset was the primary software used for simulations. It is 
generally assumed that detailed modeling of the buildings and it environment are required 
for an accurate prediction of the signal contours since the signal contour prediction by 
Urbana will only be as good as the accuracy of the information provided. Urbana does 
not take into account the mobility of objects in the models and natural attenuation due to 
atmospheric conditions like rain and fog. After the simulations were completed, Matlab 
and Microsoft Excel were used to analyze data and display tabular results. However, the 
field strength contour plots were generated by Urbana. 

Typical cell phone frequencies of 850 MHz and 900 MHz were selected for this 
research. Several other input parameters were also varied in the process of simulation, 
such as, building material, observation areas, transmitting power, and antenna patterns, to 
obtain sufficient data to satisfy the objectives for this research. 
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B. CONCLUSIONS 


The results and data analysis from the simulation have shown that varying 
frequeney over the speeified range did not show signifieant impaet on the overall results. 
However, the building material ean eontribute some impaet on the signal propagating in 
the environment. The average power for buildings with PEC material is eonsistently 
higher by 9 dBm in eomparing to the ease of buildings with eonerete material. This was 
primarily due to the faet that the loeation of the observation grid happened to be in an 
area of strong refleetions. Materials that have higher transmission ooeffieients result in 
more diffuse signal distribution. 

Based on the simulation results and data shown in Table 4, we ean eonelude that 
adding windows will not signifieantly ehange the signal distribution on average. It is 
important to emphasize that this statement applies to average values over the observation 
grid. There are some points that inerease or deerease substantially, but taken on the 
average, the levels remain about the same. Therefore it may not be neeessary to build a 
detailed model with the eorreet building material in order to obtain good results from the 
simulations. 

The seeond objeetive was to determine whieh antenna options would provide the 
maximum eoverage over a speeified seetor. The results from these simulations and data 
shown in Table 7 indieate that it is no more advantageous to utilize two antennas with 
lower transmitting power and higher antenna gain than a single wider beamwidth antenna 
that uses a higher transmitting power. 

C. FUTURE WORK 

A more detailed profiling of an urban environment sueh as in downtown Los 
Angeles would be useful. Model details should inelude billboards, lamp posts, street 
signs, trees, and parked ears to fully assess their impaet on simulation results. 

Future work ean also inelude validation and verifieation of the simulation results 
through aetual measurements. 
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Multiple diffraction should be investigated. Unfortunately Urbana only does 
“over the roof’ diffraction (not multiple diffraction from vertical edges). Reference [4] 
indicates that multiple diffraction from vertical edges can be significant. 
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